arrows(x0=5,y0=8,x1=5,y=2, lwd=2)
arrows(x0=2,y0=5,x1=8,y=5, lwd=2)
arrows(x0=8,y0=5,x1=2,y=5, lwd=2)
text(0.75,5,"Organization of \n Production", font=2)
text(5,8.5,"Trade Policy", font=2)
par(mar=c(0,0,0,0))
x <- c(-1,-2,-3)
y <- c(-1,-2,-3)
plot(x~ y, xlim=c(0,10), ylim=c(0,10), xaxt="n", yaxt="n",
xlab="", ylab="", bty="n")
arrows(x0=5,y0=2,x1=5,y=8, lwd=2)
arrows(x0=5,y0=8,x1=5,y=2, lwd=2)
arrows(x0=2,y0=5,x1=8,y=5, lwd=2)
arrows(x0=8,y0=5,x1=2,y=5, lwd=2)
text(0.75,5,"Organization of \n Production", font=2)
text(5,8.75,"Trade Policy", font=2)
text(5,8.25,"Protectionist")
text(5,1.25,"Liberal")
text(5,2,"Liberal")
text(5,1.75,"Liberal")
par(mar=c(0,0,0,0))
x <- c(-1,-2,-3)
y <- c(-1,-2,-3)
plot(x~ y, xlim=c(0,10), ylim=c(0,10), xaxt="n", yaxt="n",
xlab="", ylab="", bty="n")
arrows(x0=5,y0=2,x1=5,y=8, lwd=2)
arrows(x0=5,y0=8,x1=5,y=2, lwd=2)
arrows(x0=2,y0=5,x1=8,y=5, lwd=2)
arrows(x0=8,y0=5,x1=2,y=5, lwd=2)
text(0.75,5,"Organization of \n Production", font=2)
text(5,8.75,"Trade Policy", font=2)
text(5,8.25,"Protectionist")
text(5,1.75,"Liberal")
text(1.5,5,"State")
text(1.5,5,"Market")
text(1.5,4,"State")
text(1.5,4.5,"State")
text(2),4.5,"State")
text(2,4.5,"State")
text(3,4.5,"State")
text(2.5,4.5,"State")
par(mar=c(0,0,0,0))
x <- c(-1,-2,-3)
y <- c(-1,-2,-3)
plot(x~ y, xlim=c(0,10), ylim=c(0,10), xaxt="n", yaxt="n",
xlab="", ylab="", bty="n")
arrows(x0=5,y0=2,x1=5,y=8, lwd=2)
arrows(x0=5,y0=8,x1=5,y=2, lwd=2)
arrows(x0=2,y0=5,x1=8,y=5, lwd=2)
arrows(x0=8,y0=5,x1=2,y=5, lwd=2)
text(0.75,5,"Organization of \n Production", font=2)
text(5,8.75,"Trade Policy", font=2)
text(5,8.25,"Protectionist")
text(5,1.75,"Liberal")
text(2.25,4.5,"State")
text(7.75,4.5,"Market")
x <- (6, 8, 10, 12, 14)
x <- c(6, 8, 10, 12, 14)
y <- x
plot(y ~ x)
plot(y ~ x, lty=1)
plot(y ~ x, lty="line")
plot(y ~ x, lty=2)
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13))
plot(y ~ x, lty=2, xlim=c(6,12), ylim=c(6,12))
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
lty=1)
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type=1)
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type=1)
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
lty="l")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l")
x <- c(6, 8, 10, 12, 14)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=2)
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="line", lwd=2)
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=2)
,
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=2)
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=2, ylab = "Net worth in 2005 (in logarithm)")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=2, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=3, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=3, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
x <- c(6:14, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=3, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
x <- c(6:14, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=3, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=3, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
?segments
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="lty", lwd=3, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
?segments
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
lty="lty", lwd=3, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=3, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
segments(x0=6, y0=6, x1 = 14, y1 = 14)
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
segments(x0=6, y0=6, x1 = 14, y1 = 14, lwd=5)
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)")
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
?axis
axis(side=1, at=c(6,8,10,12))
axis(side=2, at=c(6,8,10,12))
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main="Figure 1: Wealth in 1999 and wealth in 2005, Thailand")
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"))
x <- c(6:14, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"))
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(6,8,10,12))
axis(side=2, at=c(6,8,10,12))
x <- c(6:14, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"))
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(6,8,10,12))
axis(side=2, at=c(6,8,10,12))
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"))
x <- c(6:14, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"))
plot(lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"))
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"))
x <- c(6:14, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"),
axes=F)
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(6,8,10,12))
axis(side=2, at=c(6,8,10,12))
x <- c(6:13, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"),
axes=F)
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(6,8,10,12))
axis(side=2, at=c(6,8,10,12))
axis(side=1, at=c(5,6,8,10,12))
x <- c(6:13, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"),
axes=F)
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(5,6,8,10,12))
axis(side=2, at=c(5,6,8,10,12))
axis(side=1, at=c(5,6,8,10,12,14))
axis(side=2, at=c(5,6,8,10,12,14))
x <- c(6:13, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"),
axes=F)
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(5,6,8,10,12,14))
axis(side=2, at=c(5,6,8,10,12,14))
mtext(10,10,"P")
text(10,10,"P")
?mtext
?text
text(x=10,y=10,"P")
text(x=7.5,y=7.5,"P")
text(x=7.5,y=7,"P")
x <- c(6:13, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"),
axes=F)
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(5,6,8,10,12,14))
axis(side=2, at=c(5,6,8,10,12,14))
text(x=7.5,y=7,"P",cex=2)
x <- c(6:13, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"),
axes=F)
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(5,6,8,10,12,14))
axis(side=2, at=c(5,6,8,10,12,14))
text(x=7.5,y=7,"P",cex=1.75)
x <- c(6:13, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"),
axes=F)
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(5,6,8,10,12,14))
axis(side=2, at=c(5,6,8,10,12,14))
text(x=7.5,y=7,"P",cex=1.75)
x <- c(6:13, by=0.1)
y <- x
plot(y ~ x, lty=2, xlim=c(6,13), ylim=c(6,13), xaxt="n", yaxt="n",
type="l", lwd=1, ylab = "Net worth in 2005 (in logarithm)",
xlab = "Net worth in 1999 (in logarithm)",
main=expression("Figure 1: Wealth in 1999 and wealth in 2005, Thailand"),
axes=F)
segments(x0=5, y0=5, x1 = 14, y1 = 14, lwd=3)
axis(side=1, at=c(5,6,8,10,12,14))
axis(side=2, at=c(5,6,8,10,12,14))
text(x=7.5,y=7,"P",cex=1.75)
x <- seq(1,48,929,184,2929)
x <- c(1,48,929,184,2929)
y <- c(29,2,28,10,5)
cor(x,Y)
cor(x,y)
summary(lm(x ~ y))
x <- seq(0,100, by=0.1)
y <- x^(10)
plot(y ~ x)
x <- seq(0,100, by=0.1)
y <- x^(20)
plot(y ~ x)
x <- seq(0,100, by=0.1)
y <- x^(50)
plot(y ~ x)
plot(y ~ x, lty=1)
plot(y ~ x, lty=2)
x <- seq(0,100, by=0.1)
y <- x^(50)
plot(y ~ x, lty=2)
plot(y ~ x, lty="l")
plot(y ~ x, lty="line")
x <- seq(0,100, by=0.001)
y <- x^(50)
plot(y ~ x)
x <- seq(0,100, by=0.001)
y <- x^(50)
plot(y ~ x)
x <- seq(0,100, by=0.001)
y <- x^(50)
plot(y ~ x, lty="solid")
x <- seq(0,100, by=0.1)
y <- x^(50)
plot(y ~ x, lty="solid")
x <- seq(0,100, by=0.1)
y <- x^(50)
plot(y ~ x, type="l")
x <- seq(0,100, by=0.1)
y <- x^(100)
plot(y ~ x, type="l")
x <- seq(0,100, by=0.1)
y <- x^(100)
plot(y ~ x, type="l", ylab="% of people at social function who are white")
plot(y ~ x, type="l", ylab="% of people at social function who are white")
x <- seq(0,100, by=0.1)
y <- x^(100)
plot(y ~ x, type="l", ylab="% of people at social function who are white")
plot(y ~ x, type="l", xlab="% of people at social function who are white")
plot(y ~ x, type="l", xlab="% of people at social function who are white",
ylab="probability somebody wears bf")
plot(y ~ x, type="l", xlab="% of people at social function who are white",
ylab="probability somebody wears bf", yaxt="n")
?axis
axis(2)
axis(2, lables=c("1","2","3","4","5", "6"))
?axis
?axis
axis(2, labels=c("1","2","3","4","5", "6"))
axis(2, labels=c("1","2","3","4","5", "6"), at=NULL)
x <- seq(0,100, by=0.1)
y <- x^(100)
plot(y ~ x, type="l", xlab="% of people at social function who are white",
ylab="probability somebody wears bf", yaxt="n")
plot(y ~ x, type="l", xlab="% of people at social function who are white",
ylab="probability somebody wears bf")
axis(2, labels=c("0%", "100%"), at=c(0,1e+200))
x <- seq(0,100, by=0.1)
y <- x^(100)
plot(y ~ x, type="l", xlab="% of people at social function who are white",
ylab="probability somebody wears bf")
axis(2, labels=c("0%", "100%"), at=c(0,1e+200))
plot(y ~ x, type="l", xlab="% of people at social function who are white",
ylab="probability somebody wears bf", yaxt="n")
axis(2, labels=c("0%", "100%"), at=c(0,1e+200))
plot(y ~ x, type="l", xlab="% of people at social function who are white",
ylab="probability somebody wears bf")
axis(2, labels=c("0%", "20%", "40%", "60%", "80%", "100%"), at=c(0, 2e+199, 4e+199,
6e+199, 8e+199, 1e+200))
x <- seq(0,100, by=0.1)
y <- x^(100)
plot(y ~ x, type="l", xlab="% of people at social function who are white",
ylab="probability somebody wears bf", yaxt="n")
axis(2, labels=c("0%", "20%", "40%", "60%", "80%", "100%"), at=c(0, 2e+199, 4e+199,
6e+199, 8e+199, 1e+200))
466_168
466+168
1381+125
1500+650+450
415.21
3544.534
3544.534/415.21
415*.6
3544.534/250
16.5/15
1.1*8.5
30046/(11718 + 4707)
3.147774947/.415
7.585/1.829285
404647020/3147774947
250/315
11718+4707
16425/12649
800/415
1.298522*1.927711
404647020/30046
415000000*.4
1.66e+08
166000000/(11718+4707)
65-18
1500000/47
1000000/47
20*(1.03)
20.6*1.03
(20*(1.03))^2
20*(1.03^2)
20*(1.03^300)
# 6 million 300 years later
100*(1.03^300)
100*(1.05^300)
100*(1.04^300)
100*(1.035^300)
100*(1.036^300)
100*(1.037^300)
100*(1.04^300)
20*(1.04^300)
20*(1.05^300)
20*(1.045^300)
100*(1.045^300)
100*(1.04^300)
100*(1.03^300)
100*(1.035^300)
20*(1.035^300)
20*(1.04^300)
20*(1.54^300)
20*(1.05^300)
20*(1.045^300)
24*(1.045^300)
24*(1.04^300)
24*(1.045^300)
24*(1.0425^300)
24*1.0425
24*1.0425
25.02*1.0425
24*(1.04^30)
24*(1.05^30)
24*(1.40^30)
100*(1.40^30)
100*(1.50^30)
100*(1.45^30)
100*1.45
145*1.45
6*(1.03*189)
6*(1.10*189)
6*(1.100*189)
6*(1.10000*189)
6*(1.5*189)
6*(10*189)
6*(10*189)
6*10
6*(1000*189)
6*(10000*189)
6*(1.5*18.9)
6*(2*18.9)
6*(10*18.9)
6*(1000*18.9)
6*(10000*18.9)
6*(15000*18.9)
getwd()
setwd("/Users/dinohadzic/Desktop")
data <- read.csv("Votes.csv")
names(data)
dim(data)
m1 <- glm(Vote_Yes ~ Party_D, data=data, family = binomial(link = "logit"))
summary(m1)
names(data)
mRegion <- glm(Vote_Yes ~ Region_N, data=data, family = binomial(link = "logit"))
mParty <- glm(Vote_Yes ~ Party_D, data=data, family = binomial(link = "logit"))
summary(mParty)
mRegion <- glm(Vote_Yes ~ Region_N, data=data, family = binomial(link = "logit"))
summary(mRegion)
mPartyRegion <- glm(Vote_Yes ~ Party_D + Region_N, data=data, family = binomial(link = "logit"))
summary(mPartyRegion)
setwd("/Users/dinohadzic/Desktop/Replication/OA.D")
setwd("/Users/dinohadzic/Desktop/Replication/OA.D")
data <- read.csv("OA_D_data.csv")
setwd("/Users/dinohadzic/Desktop/Replication/OA.D")
data <- read.csv("OA_D_data.csv")
names(data)
# Table OA.D1: Explaining Voting Transparency
# Model 1
mod1A <- glm(RCV_SOP ~ Average_ballot_under_control_of_party + Size_of_the_chamber +
Single_member_district_seats + Legislature_can_remove_executive + Bicameralism +
Directly_elected_members, data=data, family = binomial(link = "logit"))
# Table OA.D1: Explaining Voting Transparency
# Model 1
mod1 <- glm(RCV_SOP ~ Average_ballot_under_control_of_party + Size_of_the_chamber +
Single_member_district_seats + Legislature_can_remove_executive + Bicameralism +
Directly_elected_members, data=data, family = binomial(link = "logit"))
summary(mod1)
nobs(mod1) # 85
AIC(mod1) # 109.732
BIC(mod1) # 126.8306
logLik(mod1) # -47.86602
# Model 2
mod2 <- glm(RCV_SOP ~ Average_ballot_under_control_of_party*Size_of_the_chamber +
Single_member_district_seats + Legislature_can_remove_executive + Bicameralism +
Directly_elected_members, data=data, family = binomial(link = "logit"))
summary(mod2)
nobs(mod2) # 85
AIC(mod2) # 109.0411
BIC(mod2) # 128.5823
logLik(mod2) # -46.52055
names(data)
# Model 3
mod3A <- glm(RCV_SOP ~ Average_ballot_under_control_of_party + Size_of_the_chamber +
Single_member_district_seats + Legislature_can_remove_executive + Bicameralism +
Directly_elected_members + Polity_IV_index, data=data, family = binomial(link = "logit"))
# Model 3
mod3 <- glm(RCV_SOP ~ Average_ballot_under_control_of_party + Size_of_the_chamber +
Single_member_district_seats + Legislature_can_remove_executive + Bicameralism +
Directly_elected_members + Polity_IV_index, data=data, family = binomial(link = "logit"))
summary(mod3)
nobs(mod3) # 85
AIC(mod3) # 108.2473
BIC(mod3) # 127.7885
logLik(mod3) # -46.12363
# Model 4
mod4 <- glm(RCV_SOP ~ Average_ballot_under_control_of_party*Size_of_the_chamber +
Single_member_district_seats + Legislature_can_remove_executive + Bicameralism +
Directly_elected_members + Polity_IV_index, data=data, family = binomial(link = "logit"))
summary(mod4)
nobs(mod4) # 85
AIC(mod4) # 107.8464
BIC(mod4) # 129.8302
logLik(mod4) # -44.92318
